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Abstract 

Methanol is an amphiphilic solute whose aqueous solutions exhibit distinctive physical properties. 
The volume change upon mixing, for example, is negative across the entire composition range, 
indicating strong association. We explore the corresponding behavior of a Jagla solvent, which 
has been previously shown to exhibit many of the anomalous properties of water. We consider 
two models of an amphiphilic solute: (i) a "dimer" model, which consists of one hydrophobic hard 
sphere linked to a Jagla particle with a permanent bond, and (ii) a "monomer" model, which is a 
limiting case of the dimer, formed by concentrically overlapping a hard sphere and a Jagla particle. 
Using discrete molecular dynamics, we calculate the thermodynamic properties of the resulting 
solutions. We systematically vary the set of parameters of the dimer and monomer models and 
find that one can readily reproduce the experimental behavior of the excess volume of the methanol- 
water system as a function of methanol volume fraction. We compare the pressure and temperature 
dependence of the excess volume and the excess enthalpy of both models with experimental data 
on methanol-water solutions and find qualitative agreement in most cases. We also investigate the 
solute effect on the temperature of maximum density and find that the effect of concentration is 
orders of magnitude stronger than measured experimentally. 
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I. INTRODUCTION 



Aqueous solutions of alcohols have attracted both experimental and theoretical attention 
on account of their ubiquity and importance in the medical, personal care, transportation 
(e.g., antifreeze, fuels), and food industries, among others Methanol is a simple 

example of an amphiphilic organic solute, and its aqueous solutions exhibit many interesting 
nonidealities. Attaining a good understanding of this simple case is therefore a natural 
starting point for studies of more complex solutes in water, such as higher alcohols or 
proteins. 

Because of the increasing availability of expanded computing power, simulations have 



become an important research tool in studying aqueous methanol solutions 
optimized potential for liquid simulation (OPLS) 



12| is frequently used to mode 
molecules. To represent water as a solvent, the SPC/E 13|, TIP3P [l^, TIP4P 



17H22|. The 
alcohol 



T1P5P 



15| . and 



16| models are frequently used. On the other hand, many thermodynamic properties 



of water can be reproduced using soft-core spherically symmetric potentials, one of the most 



important of which is the Jagla model 



23|. The Jagla potential has a hard core and a linear 



repulsive ramp, and contains two characteristic length scales: a hard core a and a soft core h. 
For a range of parameters, the Jagla model exhibits a water-like 2J] cascade of structural. 



transport, and thermodynamic anomalies 
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29j. Buldyrev et al. in 2007 |30j found 



that the Jagla solvent exhibits key water-like characteristics with respect to hydrophobic 
hydration, suggesting that the water-like characteristics of the Jagla solvent extend beyond 
the pure fluid. Here, we explore this analogy further, considering the properties of solutions 
of amphipathic solutes. 

We focus on the properties of the excess volume and the excess enthalpy, and on how 
amphiphilic solutes affect the temperature of maximum density (Tmd) of a solution. At 
T = 298K and P = O.lMPa the excess volume and the excess enthalpy are negative across the 
entire range of methanol concentrations (both quantities will be defined rigorously below). 
The strongest effect occurs at a methanol volume fraction of <Apax = 59.7%, at which the 
negative excess volume deviates from additivity by —3.57% |5|]. As T increases and P 
decreases, the excess volume becomes more negative and the extremal point shifts (see 
Fig. [Hand Fig. [2]). Most solutes tend to suppress water's non-idealities, and hence they 
lower the Tmd- Some amphiphilic solutes behave differently: ethanol and t-butanol have a 
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marked non-monotonic effect (whereby the Tmd of the solution first increases with respect 
to that of water upon solute addition, but decreases for more concentrated solutions), and 
methanol has a very mild non-monotonic effect 31|. 

As a first step , in this paper, our goal is to model a simple amphiphilic solution that 
mimics the properties of the methanol-water system. We use a hard sphere with parameter- 
ized diameter au to model the hydrophobic methyl group, and a Jagla particle to model the 
hydrophilic hydroxyl group. We link the hard sphere and the Jagla particle using a bond of 
adjustable length to model the amphiphilic solutes. 

In Section II of this paper we describe in detail our models and simulation methods. 
In Section III we list and analyze the simulation results, and in its four subsections we 
discuss the different parameter effects, the temperature and pressure dependence of the 
excess volume, the behavior of the excess enthalpy, and how the solute concentration affects 
the temperature of maximum density. In Section IV we list our main conclusions. 



II. MODEL AND METHODS 



A. Dimer model for amphiphilic solutes 

To model the simple amphiphilic solute methanol, we separate CH3OH into the methyl 
(CH3) and the hydroxyl (OH) groups. We model CH3 as a hard sphere and OH as a Jagla 
particle. There is one bond between the hard sphere and the Jagla particle. We call this 
model the dimer model. To model the solvent, H2O, we use the Jagla particle. The following 
interactions are included: there is a Jagla potential between Jagla solvent particles, a Jagla 
potential between the Jagla solvent and a dimer's Jagla particle, and a Jagla potential 
between two dimer's Jagla particles, all of which are denoted by ?7jj(r) [Fig. |3]^a)]. The 
interaction between the hard spheres is modeled by a hard-core potential U}ii{{r), and the 
interaction between the hard spheres and the Jagla particles is modeled by a hard core 
potential Uj}i{r). We model the covalent bond with a narrow square well potential bounded 
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by two hard walls Uhondif)- The interaction potentials are 



-Uo 



oo r < a 

(Uo+UR)(b-r) ^ ^ , 

^ , — - a < r < 

b—a 

b < r < c 

" c—u 

r > c 



" c—b 



where a is the hard core diameter, h = 1.72a is the soft core diameter, c = 3a is the range of 
attractive potential, Uo is the maximum attractive energy and f/^j = 3.56f/o is the maximum 
repulsive energy, and 

oo r < au 
r > au 

where au is the diameter of the hard sphere, and 



UHHir) 



(2) 



UjHir) 



oo r < ajM 
r > ajM 



where aju 



and 



oo r < C - I 

c-f<r<C+ 

oo r > C + I 



(3) 



(4) 



The hard core diameter au and the average length of the covalent bond ( are used as 
adjustable parameters in order to achieve agreement between the excess volume of the model 
solution and the experimental results of methanol-water solutions at ambient conditions. In 
all our simulations, we use the same set of Jagla potential parameters, b = 1.72a, c = 3a, 
and Ur = 3.56Uo [29|. We use reduced units in terms of length a, energy Uo, and particle 
mass m. For temperature we use units of Uo/kB] for pressure we use units of Uo/a^; and for 
volume we use units of a^. 



B. Monomer model for amphiphilic solutes 

The best agreement between the dimer model and the experimental data occurs when 
bonds are short (for a more detailed discussion, see Sec. II), i.e., the hard sphere that models 
the methyl group and the Jagla particle that models the hydroxyl group almost overlap. For 
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this reason we also consider a methanol model in which the overlap is complete, and the 
bond length vanishes. This leads to a spherically symmetric potential that superimposes the 
Jagla particle and the hard sphere. In this "monomer" model we introduce the interaction 
potentials between "methyl" monomers, UMui^) [Fig- E^b)], between monomer and Jagla 
particle, UjM{r) [Fig. [31(c)], and between Jagla particles t/jj(r) [Fig. [31^ a)]. In this case the 
interaction formulae are 

,00 r < ttM 
UMuir) ={ , (5) 

Ujjir) r > au 



and 



,00 r < ajM 
Umir) ={ , (6) 

Ujj{r) r > ajM 



where Ujj{r) is defined by Eq. (1). 



C. Simulation details and analysis methods 

For our simulations we use the discrete molecular dynamics (DMD) algorithm. With 



DMD we approximate a continuous potential 



steps. We use the same scheme as in Ref. 



Dy a discrete potential made up of a series of 



29| . Our simulation consists of a fixed number 



= 2000 particles in a cubic box with periodic boundaries. We denote the solute mole 
fraction by x. Since the dimer contains two particles and the monomer one, the number of 
solute molecules Ns{x) is 

N.(.) - ^ (7) 

in the dimer system, and 

Ns{x) = Nx (8) 
in the monomer system. The number of the solvent particles Nj{x) is 

NAx) = ^'^-"> (9) 

' x + 1 ^ ' 

in the dimer system, and 

Nj{x) = N{l-x) (10) 
in the monomer system. The total number of molecules in the system Nt{x) is 



in the dimer model, and 



Nt{x) = N (12) 

in the monomer modeL The volume occupied by Nj{x) pure Jagla solvent particles before 
mixing, Vj(x), and the volume occupied by Ns{x) pure solute molecules before mixing, Vs{x), 
at the given temperature and pressure, are given by 

Nj(x) 

Vj{x) = ^K„ix(0) (13) 

and 

Vsix) = ^K.ix(l), (14) 

where Vmix(x) is the volume of the mixture with mole fraction x. 

We define the excess volume of the solution with respect to the ideal mixture as 

AV= , - 1 (15) 

Vj{x) + Vsix) ^ ^ 

If the excess volume AV^ is negative, the volume of the solution is less than the volume of 
the ideal mixture. If it is positive, the system expands after mixing at fixed temperature 
and pressure. In most contexts, we use the volume fraction 

rather than the mole fraction x to express different solute concentrations of solutions. 



We compare our simulation results with the data from experiments 



4, l5| , where the excess 



volume is expressed in terms of AY = ^ — "* " , where Um is the number of the moles 
of methanol, n,,, is the number of moles of water, the mole fraction is x = — — , and 
and Vm are the molar volumes of water and methanol, respectively, at specific temperature 
and pressure conditions. The conversion formulas between AV and AY, and if and x are 

Ay , , 

A ^ = r. ^ 17 

XVrr,^ + 1 - x)v,„ 



and 



^=x + (l-x)^' ^^^^ 



Density is an important system property. For this purpose, we assume that the Jagla 
particles and the solute particles correspond to the same number of water and methanol 
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molecules in a pure solution, respectively, and express the density of the pure solute in 
terms of the density ratio 

32 

P=% (19) 

where vj = -^j^ and Vg = w^i) volume per particle of the pure solvent and the pure 

solute, respectively. We compare the simulation with the experimental number p = 0.79. 
The excess enthalpy is usually defined as 

where is the total enthalpy of moles of pure methanol, and is the total enthalpy 
of moles of pure water under specific temperature and pressure conditions. To put the 
enthalpy comparison on the same footing as the excess volume data, we choose to also report 
the excess enthalpy on a volumetric basis and define the excess enthalpy per volume as 

Hi.) - ^Hjl) - f^g(O) 

^"■^ VOW + • 

where H{x) is the enthalpy of the system with a mole fraction x. The conversion formula is 

AHs = - . (22) 

In our simulation, we measure AHg in units of Uo/a^- In order to compare our results 
with experimental data, we need to convert our units into J/cm^ = MPa. In accordance 

n 

with Ref. [27|, we use Uo = 4.75KJ/mol and a = 2.7 x 10 ^cm. Then we convert by simply 
multiplying our simulation results by 4.008 x 10^. 



III. RESULTS AND DISCUSSION 



A. Effects of the parameters on model behavior 

Because there are several parameters in our dimer model, we first investigate how these af- 
fect the simulations, searching for a set of parameters that can best model ambient methanol. 
Since the goal is to explore the excess properties of the model, vis a vis real methanol-water 

n 

solution, we compare our simulation results with the results reported in Ref. [5| concerning 
the excess volume at T = 298K and P = 0.1 MPa. We set our simulation temperature and 
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pressure at T = 0.5 and P = 0.02, and we change the diameter of the hard spheres that 
model the methyl group in methanol. In Fig. HI we show the excess volume of the solutions 
with different hard core diameters of the hard spheres om, for a fixed bond length ( = 0.9a, 
across the entire range of amphiphilic solute volume fraction. If ayi < 1.3a, we cannot 
reproduce the volume reduction over the entire range of volume fraction. If om > 1.3a, the 
mixture shrinks over the entire range of volume fractions. When the hard core diameter 
increases, the volume fraction corresponding to maximum volume reduction (henceforth re- 
ferred to as maximum reduction volume fraction) decreases. However, for this value of ( we 
cannot reproduce the experimentally-observed maximum reduction volume fraction. 

We next show, in Fig. El how bond length affects the excess volume. We choose three 
different hard sphere diameters and vary the bond length. For cm = l-5a [Fig. El^a)], as the 
bond length increases, the excess volume becomes more negative and the maximum reduction 
volume fraction increases, but its value is always larger than in experiments. For om = 1.6a 
[Fig. El^b)] and om = l-7a [Fig. EJ^c)], the excess volume and the maximum reduction volume 
fraction exhibit the same trend as the bond changes as seen for au = 1.5a. However, their 
maximum reduction volume fraction is smaller and closer to the experimental data. Thus, 
we can see that the maximum reduction volume fraction (/?max and the corresponding AV^ax 
agree well with the experimental observations for selected bond lengths; i.e., for au = l-6a 
and C = 0.35a, v^max = 63.6%, A\4iax = —3.67%, and for an = l-7a and ( = 0.125a, 
</2max = 58.717%, A\4iax = —3.55%. In addition, over the entire range of volume fractions 
the excess volume agrees quantitatively with experiment for these two sets of parameters. 

In the dimer model, we achieve the closest agreement with the experimental data at 
qm = 1-7 a and ( = 0.125a. At this set of parameters, the hard sphere and the Jagla particle 
nearly overlap. The dimer model can thus be modeled as a monomer with a hard core 
au ~ 1.7a, as shown in Eqs. (5) and (6). In order to achieve the closest agreement with 
the experimental data, we vary the single parameter, au in this model (Fig. [6]). For all 
cases of au, the reduction property can be reproduced across the entire range of volume 
fractions. At au = 1.73a and au = 1.75a, the trend of the excess volume is very similar to 
the experimental results. For au = 1. 74a, the maximum reduction volume fraction is close 
to the experimental data, and in the range of volume fractions smaller than the maximum 
reduction volume fraction, the excess volume is also quantitatively reproduced. Overall, the 
monomer model can reproduce the excess volume curve qualitatively but is not as accurate 
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as the dimer model. 

At T = 298K and P = 0.1 MPa, the methanol density is 0.78663 g/cc [sl, and its ratio 
with respect to water is 0.79. We have explored the relative density of the neat amphiphilic 
solutes with respect to a pure Jagla solvent when we change the parameters. The density 
increases when we shorten the bond length, when the diameter of the hard sphere om in 
the dimer model decreases, and when the hard core of the monomer om decreases. From 
the above- discussed results we know that when au = ^-^cl and ( = 0.125a, we achieve the 
best agreement between our simulation results and the experimental results. However, the 
density ratio of the pure solute for this set of parameters is 0.66, less than the density ratio of 
real methanol at T = 298 K and P = 0.1 MPa. In the monomer model, when the hard core 
diameter of monomer is = 1.74a, the density ratio is 0.64, also less than the experimental 
value. 

B. Temperature and pressure dependence of the excess volume 

Experimentally, when the temperature change is small, e.g., 5K, the change of the max- 
imum reduction volume fraction is negligible and there is a small increase in AV as the 
temperature increases [Fig. [T](a)]. Over large temperature intervals, e.g., 50K, the excess 
volume becomes increasingly negative and the maximum reduction volume fraction increases 
slightly as the temperature increases [Figs. [l](b) and[T|^c)]. We use our two models to check 
the temperature and pressure dependence of the excess volume. Choosing the parameters 
that can best reproduce the experimental results of the excess volume at T = 298 K and 
P = 0.1 MPa, we set au = 1.7a and ( = 0.125a for the dimer model. 

We calculate the temperature dependence of the excess volume at pressures P = 0.02 
and P = 0.1 in the dimer model [Fig. [7]. At P = 0.02 1.0 atm), we first measure the 
excess volume at a series of temperatures separated by 0.01 in the range of T = 0.5-0.6 
[Fig. El^a)]. The excess volume becomes slightly more negative and the change of the maxi- 
mum reduction volume fraction is negligible as the temperature increases. When we enlarge 
the temperature interval to 0.05 [Fig. EJ^b)], the excess volume becomes increasingly nega- 
tive and the maximum reduction volume fraction increases noticeably as the temperature 
increases. At P = 0.1 [Fig. [Tl^c)], we find approximately the same results. However, when 
the temperature is approximately 0.5, the excess volume of the solution with less than 15 
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percent of amphiphilic solutes goes to zero or becomes positive, indicating that in our model 
the volume does not decrease, and can even expand at these conditions. We will discuss 
this "bump" in the excess volume curve later. The temperature dependence at pressures 
P = 0.02 and P = 0.1 using the Cm = 1.74a monomer model (not shown) are comparable 
to those of the dimer model. 

We next compare the calculated pressure dependence of the excess volume with the 
experimental results. In Fig. [2l we report the experimental pressure dependence of the excess 
volume at two different temperatures, T = 323. 15K and T = 283. 15K, using two different 
pressure intervals. The results at T = 323. 15K cover pressures from P = 0.1 MPa to P = 
13.5MPa, and those at T = 283. 15K cover pressures from P = O.lMPa to P = 205MPa. We 
can see that, as the pressure increases, the excess volume becomes less negative. Regarding 
the maximum reduction volume fraction of the excess volume, in Fig. [2|^a), its value does 
not noticeably change with temperature, however in Fig. |2]^b), we can clearly see that as the 
pressure increases, the maximum reduction volume fraction does increase. 

In our simulation, we fix the temperature at T = 0.5 and calculate the excess volume 
at pressures from P = 0.02 to P = 0.1. Our simulation results for the dimer model (see 
Fig. E]) agree with the experimental results quite well, namely, as the pressure increases the 
excess volume becomes less negative and the maximum reduction volume fraction increases. 
Regarding the positive "bump" at the small volume fraction of amphiphilic solutes, we 
are not aware of experimental data for this range of volume fractions, so this result is 
a prediction, namely that at high pressures and low temperatures, dilute methanol- water 
solutions have a positive excess volume. The simulation results for the best monomer model 
are quite similar to those for the dimer model except that the positive "bump" is smaller. 

C. Excess enthalpy 

The excess enthalpy in methanol-water solutions has been measured at T = 298 K, and 
P = 0.1 MPa jj]. The excess enthalpy is negative (exothermic mixing), consistent with 
the picture of strong association that follows from the volumetric behavior. The maximum 
reduction occurs at a volume fraction v^max = 42%, differing from that of the excess volume 
at these specific conditions. At pressures P = O.lMPa, P = 20MPa, and P = 39MPa, the 
excess enthalpy becomes less negative as the temperature increases. The maximum reduction 
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volume fraction increases as the temperature increases. In some temperature ranges, e.g., 
T = 278.15K-T = 298. 15K, the maximum reduction volume fraction actually increases as 
)ressure increases, and the excess enthalpy becomes more negative as the pressure increases 

ml 

In the dimer model, when om > 1.0a, all the cases can reproduce qualitatively the 
enthalpy of mixing rather well. As au increases, the excess enthalpy becomes more negative 
and the maximum reduction volume fraction becomes smaller, following the same trend as 
in experiments. The magnitude of the excess enthalpy for om = l-7a and au = 1.8a can be 
made close to the experimental results, but for the best model developed in Section III. A, 
the agreement is only qualitative. 

We report the bond length dependence of the excess enthalpy at a^s = l-5a [Fig. IHI^a)], 
au = l-6a [Fig. [Uj^b)], and au = l-7a [Fig. ISt^c)]. For au = l-5a and ans = 1.6a, as 
the bond length increases, the excess enthalpy becomes more negative. For om = 1.7a, 
the excess enthalpy first becomes more negative and then less negative as the bond length 
increases. For all three values of au, the maximum reduction volume faction becomes larger 
as the bond length increases and the maximum enthalpy reduction is smaller than the 
experimental value. For the monomer model, when au > 1.6a, the excess enthalpy can also 
be qualitatively reproduced, and as au increases it becomes more negative and the maximum 
reduction volume fraction increases. The value of the excess enthalpy from simulations is 
smaller than the experimental values, but the magnitude is comparable. 

We return to the dimer with au = l-7a and ( = 0.125a and examine the temperature de- 
pendence. At P = 0.02 [Fig. [TOT a)]. our simulation results contradict the experimental data: 
the excess enthalpy becomes increasingly negative as the temperature increases. However, 
the increasing maximum reduction volume fraction with increasing temperature does agree 
with the experimental trend. At P = 0.1 [Fig. [TOT b)]. the excess enthalpy changes only very 
slightly with temperature, although at the highest T the magnitude of the excess enthalpy 
decreases slightly with T. If we increase the pressure to P = 0.15 [Fig. fTOT c)]. we can see 
that as the temperature increases, the magnitude of excess enthalpy decreases with T and 
the maximum reduction volume fraction increases, in agreement with experimental trends. 

In our simulation of the pressure dependence of the excess enthalpy at T = 0.5 (Fig. [TTj) . 
as the pressure increases, the excess enthalpy becomes more negative and the maximum 
reduction volume fraction becomes larger, which agrees well with experimental observations 
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of methanol-water solutions. 



D. Effect on the temperature of maximum density 

At sufficiently low temperatures and pressures, the density of liquid water exhibits a 
maximum with respect to temperature at fixed pressure. For example, liquid water has 
a maximum density at T = 277K and P = O.lMPa. If we add solutes to water, the 



temperature of the maximum density changes. According to Ref . [3l| , the Tmd of a methanol 
solutions reaches its maximum at around x = 0.6% mole fraction and then decreases slightly 
as the mole fraction increases, reaching 269K at x = 5%. This non-monotonic behavior has 



been explained by Chatterjee et al. using a statistical mechanical model of water [32 1. 

We explore the change of the Tmd with concentration for our model. We define the 
change of the temperature of maximum density as AT = Tmds — ^mdj where Tmds is Tmd 
of the solution and Tmdj is the Tmd of the pure Jagla liquid, at the given pressure. We have 
carried out simulations at three different pressures for the case of dimer with om = l-7a 
and ( = 0.125a. Our results are shown in Table I and Fig. [121 We find that AT is always 
negative and its absolute value increases monotonically with solute mole fraction. If we 
increase pressure for the same mole fraction, the change of Tmd decreases in magnitude, but 
it never becomes positive. Moreover, we find that the decrease of Tmd in our model is orders 
of magnitude stronger in methanol-water mixtures. 



IV. CONCLUSION 



Inspired by the distinctive properties of methanol-water solutions, we construct a dimer 
with a hard sphere and a Jagla particle to model amphiphilic solutes. We vary the hard 
core diameter and the bond length to achieve the best agreement between simulations and 
experiment in the excess volume. We find that the best agreement occurs for the dimer 
with Am = l-7a and ( = 0.125a, which suggests that the dimer model can be reduced 
to a monomer with a large hard core and an attractive potential that coincides with the 
attractive part of the Jagla potential. Regarding the temperature and pressure dependence 
of the excess volume, our results agree qualitatively with experimental data. Our model 
reproduces the excess enthalpy of the methanol solutions less accurately than the excess 
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volume. This is related to the fact that, in our simple model of amphiphilic solutes, we 
use the unchanged Jagla potential for the amphiphilic group. We speculate that a better 
agreement could be achieved if we varied the potential of the amphiphilic group. When we 
investigate the effect of the amphiphilic solute on the temperature of maximum density of a 
solution, we find that unlike in water-met hanol solution, the Tmd monotonically decreases 
with solute concentrations. Moreover, the effect of concentration in the model is orders of 
magnitude stronger than in experiments. 
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TABLE I: The temperature of maximum density, TmD) and the difference between the Tmd of the 
solution and that of the pure solvent at the given pressure, AT, for solutions with different mole 
fractions of amphiphilic solutes in Jagla solvents, at different pressures. Results are for the dimer 
solute model with au = 1.7a, ( = 0.125a. 
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FIG. 1: (color online) Experimental data P, [l^] on the temperature dependence of the excess 
volume, as a function of the volume fraction of methanol, at different pressures (Dotted lines 
are polynomial fits of experimental data): (a) P = O.lMPa, (b) P = 7MPa, and (c) 13.5MPa. 
Note that in (a) temperature interval between two measurements is 5K, while the temperature 
difference in (b) and (c) is 50K. The excess volume is the relative difference between the volume 
occupied by the mixture and the sum of the volumes of the pure components before mixing, at 
fixed temperature and pressure 
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FIG. 2: (color online) Experimental data [7|, llO|] on the pressure dependence of the excess volume 
as a function of the volume fraction of methanol, at two different temperatures (dotted lines 
are polynomial fits of experimental data): (a)T = 323. 2K and (b) T = 285. 15K. Note that the 
range of the pressure shown in (a) is smaller than that in (b). The excess volume is defined in the 
caption of Fig. [T] 
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FIG. 3: Sketch of the interaction potentials in our simulation, (a) The spherically symmetric 
"two-scale" Jagla ramp potential. The two length scales are the hard core diameter r = a, and 
the soft core diameter r = b. We study the case Ur = 3.56C/o, b = 1.72a and a long range cutoff 
c = 3a. (b) The interaction potential between two monomers, with contact distance cm- (c) The 
interaction potential between the Jagla particle and the monomer, contacting at (a + aM)/2. 
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FIG. 4: (color online) Dependence of the excess volume as a function of the solute volume fraction 
for a range of hard core diameters ayi, for the dimer model with (" = 0.9a. The simulations 
are done at T = 0.5 and P = 0.02. The experimental result (circles) are for T = 298K and 

n 

P = O.lMPa [5[, where the maximum volume fraction reduction occurs at tpmax ~ 59.7%. For the 
simulation results, when om < 1.3a, the excess volume has both positive and negative values. As 
om increases, the excess volume becomes more negative and ipmax, the maximum reduction volume 
fraction, decreases. 
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FIG. 5: (color online) Dependence of the excess volume on solute volume fraction for various solute 
bond lengths ( at three values of hard core diameter am- The simulations are done at T = 0.5 and 
P = 0.02. (a) OM = 1.5a, (b) om = 1.6a, (c) om = 1.7a. As ( increases, the excess volume, AV 
becomes more negative and the maximum reduction volume fraction (pma.x increases. Calculations 
performed using the parameter sets C = 0.35a, au = l-6a in (b) and ( = 0.125a, um = l-7a in (c) 
agree well with experiment. 
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FIG. 6: (color online) Dependence of the excess volume on solute volume fraction for various hard 
core diameter values (om), in the monomer ((" = 0) model. Simulations at T = 0.5 and P = 0.02. 
As Am increases, AV becomes more negative; the maximum reduction volume fraction '^max 
relatively insensitive to um- The excess volume at cm = 1.74a is similar to experiment, but the 
results are less accurate than for the dimer model [see Fig. [St^b) and[5|^c)]. 
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FIG. 7: (color online) Dependence of the excess volume on solute volume fraction, at various 
temperatures and three different pressures, for the dimer model with om = 1.7a, C, = 0.125a. 
(a) P = 0.02, which is comparable to atmospheric pressure. The temperature difference between 
curves is 0.01. (b) P = 0.02. The temperature difference between curves is 0.5. (c) P = 0.1. As 
the temperature increases, the excess volume becomes more negative and the maximum reduction 
volume fraction increases. In (c), when the temperature is low, there is a range of dilute mixtures 
for which l\V > 0. There are no corresponding experimental observations. 
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FIG. 8: (color online) Dependence of the excess volume on solute volume fraction for several 
different pressures, at T = 0.5, for the dimer model (om = 1.7a, C = 0.125a). As the pressure 
increases, the excess volume becomes more negative and the maximum reduction volume fraction 
increases very slightly. 



23 




) I , I I I , I , L_ 

0.2 0.4 0.6 0.8 

9 




) I , I I I , I , L_ 

0.2 0.4 0.6 0.8 

9 




0.2 0.4 0.6 0.8 

9 



FIG. 9: (color online) Dependence of the excess enthalpy on solute volume fraction for various 
values of the bond length, at three different hard core diameters cm, in the dimer model. The 
simulations are done at P = 0.02 and T = 0.05. (a) om = 1.5a, (b) om = 1.6a, (c) au = 1.7a. For 
(a) and (b), as the bond length increases, the excess enthalpy becomes more negative. For (c), the 
excess enthalpy exhibits a non-monotonic dependence on bond length. The maximum reduction 
volume fraction increases as C increases for all three cases. 
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FIG. 10: (color online) Dependence of the excess enthalpy on solute volume fraction for various 
temperatures, for the dimer model with cm = l-7a, Q = 0.125a. (a) P = 0.02. As the temperature 
increases, the excess enthalpy becomes more negative, in contrast to experiment, and the maximum 
reduction volume fraction increases ^. (b) P = 0.1. The excess enthalpy is nearly independent 
of temperature, (c) P = 0.15. When T > 0.6, the excess enthalpy decreases as the temperature 
increases and the volume fraction corresponding to maximum exothermicity shifts towards slightly 
higher values. 
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FIG. 11: (color online) Dependence of the excess enthalpy on solute volume fraction for various 
values of the pressure, at T = 0.5, for the dimer model with am = 1.7a, ( = 0.125a. As the pres- 
sure increases, the excess enthalpy increases, and the volume fraction corresponding to maximum 
exothermicity increases. 
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FIG. 12: (color online) Effect of amphipliilic solutes on Tmd of solutions. We use the dimer 
with parameters of am = 1.7a, C, = 0.125a. (a) The Tmd of the solutions at different solute 
concentrations, for P = 0.1. (b) The shift of Tmd; with respect to that of the pure solvent, AT, as 
a function of solute mole fraction, x, at three different pressures. Note that AT is always negative, 
indicating that adding solute lowers the solution's Tmd (lower temperatures are needed to observe 
negative thermal expansion). For a given mixture, the magnitude of AT decreases as the pressure 
increases. 



27 



